Fourier Analysis 



Fourier series allow you to expand a function on a finite interval as an infinite series of trigonometric 
functions. What if the interval is infinite? That's the subject of this chapter. Instead of a sum over 
frequencies, you will have an integral. 

15.1 Fourier Transform 

For the finite interval you have to specify the boundary conditions in order to determine the particular 
basis that you're going to use. On the infinite interval you don't have this large set of choices. After 
all, if the boundary is infinitely far away, how can it affect what you're doing over a finite distance? But 
see section 15.6. 

In section 5.3 you have several boundary condition listed that you can use on the differential 
equation u" = Xu and that will lead to orthogonal functions on your interval. For the purposes here 
the easiest approach is to assume periodic boundary conditions on the finite interval and then to take 
the limit as the length of the interval approaches infinity. On —L < x < +L, the conditions on the 
solutions of u" = Xu are then u{—L) = u{+L) and u'{—L) = u'{+L). The solution to this is most 
conveniently expressed as a complex exponential, Eq. (5.19) 

u{x) = e^^^, where u{-L) = e~'^^^ = u{L) = e^^^ 

This implies e^*'^^ = 1, or 2kL = 2nn, for integer n = 0, ±1, ±2, . . .. With these solutions, the other 
condition, u'{—L) = u'{+L) is already satisfied. The basis functions are then 

Un{x) = e*^"^ = e"^^^/^, for n = 0, ±1, ±2, etc. (15.1) 
On this interval you have the Fourier series expansion 

oo oo 

f{x) = 'Y^anUn{x), and {um, f) = {um-.^Yl,'^'^'^'^) ^ ^'^{'^^■''^'^) (15.2) 

— oo — oo 

In the basis of Eq. (15.1) this normalization is (um,Um) = 2L. 
Insert this into the series for /. 

oo r\ ^ oo 

f{x)= / \ Un{x) = -r {UnJ)Un{x) 

n=-oo\^^'^^/ ^^n=-oo 

Now I have to express this in terms of the explicit basis functions in order to manipulate it. When you 
use the explicit form you have to be careful not to use the same symbol (x) for two different things 
in the same expression. Inside the (un, f) there is no "x" left over — it's the dummy variable of 
integration and it is not the same x that is in the Un{x) at the end. Denote kn = Tcn/L. 

/(^) = iL ^ / dx'un{xrf{x')un{x) = ^ ^ / dx' e-^'^"^' f (x') e^''^'^ 

Now for some manipulation: As n changes by 1, changes by Akn = n/L. 

fL 



f{x) = ^ Yx c^a;'e-*^"^7(x')e*^"^ 

n=-oo -^'^ 

:L ^ e^^^"^ Afcn / da;'e-*^"^7(x') (15.3) 

" — -oo -^-L 
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For a given value of k, define the integral 



dx'e-'^''' f{x') 



If the function / vanishes sufficiently fast as x' — > oo, this integral will have a limit as L — )• oo. Call that 
limit g{k). Look back at Eq. (15.3) and you see that for large L the last factor will be approximately 
g{kji), where the approximation becomes exact as L oo. Rewrite that expression as 



fix) 



1 

27r 



E 



Akng{kr^ 



;i5.4) 



As L — 5- oo, you have Akn 0, and that turns Eq. (15.4) into an integra 

dk 



fix) 



e^^^g{k), where g{k) 



27r 



dxe *^^/(x 



:i5.5) 



The function g is called* the Fourier transform of /, and / is the inverse Fourier transform of g. 
Examples 

For an example, take the step function 




fix) 
9ik) 



1 {—a < X < a) 
(elsewhere) 

a 

dx e"'^^ 1 



then 



:i5.6) 



-ik 



hika 



2 sin ka 
k 



The first observation is of course that the dimensions check: If dx is a length then so is l/k. 
After that, there is only one parameter that you can vary, and that's a. As a increases, obviously the 
width of the function / increases, but now look at g. The first place where g{k) = is at /ca = vr. 
This value, Ti/a decreases as a increases. As / gets broader, g gets narrower (and taller). This is a 
general property of these Fourier transform pairs. 

Can you invert this Fourier transform, evaluating the integral of g to get back to /? Yes, using 
the method of contour integration this is very easy. Without contour integration it would be extremely 
difficult, and that is typically the case with these transforms; complex variable methods are essential 
to get anywhere with them. The same statement holds with many other transforms (Laplace, Radon, 
Mellin, Hilbert, etc. ) 

The inverse transform is 



dk 

277 



Akx 



2 sin ka 



dk 
2tt 



ikx ^ 



ika 



-ika 



dk 1 

2TTk 



ik 

^ik{x+a) 



Ci 



^ik{x—a) 



* Another common notation is to define g with an integral dxl\pM. That will require a cor- 
responding dkj^f^ in the inverse relation. It's more symmetric that way, but I prefer the other 
convention. 
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1. If X > +a then both x + a and x — a are positive, which implies that both exponentials vanish rapidly 
as A; — > +ioo. Push the contour C2 toward this direction and the integrand vanishes exponentially, 
making the integral zero. 

2. If —a < X < +a, then only x + a Is positive. The integral of the first term is then zero by exactly 
the preceding reasoning, but the other term has an exponential that vanishes as A; — > —zoo instead, 
implying that you must push the contour down toward —ioo. 



Jc, 27r k Jc 



C4 





^ 


C3 


( 





+i — (-l)27riRes 
27r fc=o k 



-i — ■ 27ri = 1 
27r 



C4 



The extra (—1) factor comes because the contour is clockwise. 

3. In the third domain, x < —a, both exponentials have the form e"^'^, requiring you to push the 
contour toward —ioo. The integrand now has both exponentials, so it is analytic at zero and there 
is zero residue. The integral vanishes and the whole analysis takes you back to the original function, 
Eq. (15.6). 

Another example of a Fourier transform, one that shows up often in quantum mechanics 



/(x) = e-^'/'"', so 9{k)= I rfxe-^'^^e-^'/'^' = / rfxe"^^^-^'/'"' 



The trick to doing this integral is to complete the square inside the exponent. 



-1 



1 



-ikx - xVo"^ = — [^'^ + ^^^^3; - a^k^/4 + a^k^/4] = — [{x + ika'^/2f + a^k'^/A\ 



The integral of / is now 



9{k) 



dx'e~^''^l"^ where x' = x + ika/2 

'OO 



The change of variables makes this a standard integral, Eq. (1.10), and the other factor, with the 
exponential of k'^, comes outside the integral. The result is 



g{k) = ay^e-^''''/^ 
This has the curious result that the Fourier transform of a Gaussian is* a Gaussian. 



(15.7) 



15.2 Convolution Theorem 

What is the Fourier transform of the product of two functions? It is a convolution of the individual 
transforms. What that means will come out of the computation. Take two functions /i and /2 with 
Fourier transforms gi and g2- 



dx fi{x)f2{x)e 



-ikx 



dk' 

dx I — c/i(/c')e*'^'^/2(x)e 



— ikx 



* Another function has this property: the hyperbolic secant. Look up the quantum mechanical 
harmonic oscillator solution for an infinite number of others. 
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°° dk' 



gi{k') J dx e^^'^ f2{x)e 
gi{k') j dx f2{x)e-'^^-^'^^ 
gi{k')g2{k-k') 



(15.8) 



The last expression (except for the 2n) is called the convolution of gi and g2- 

-ikx -'- 



dx fi{x)f2{x)e 



27r 



{9i*92){k) 



(15.9) 



The last line shows a common notation for the convolution o\ gi and g2- 
What is the integral of |/p over the whole line? 



dxf*{x)f{x) 



dk 

dxf*{x) I — g{k)e^'^^ 
27r 



dk 
dk 



g{k) / dxf*{x)e 



ikx 



dx f{x)e 



-ikx 



dk 

271 



g{k)g*{k) 



(15.10) 



This is Parseval's identity for Fourier transforms. There is an extension to it in problem 15.10. 
15.3 Time-Series Analysis 

Fourier analysis isn't restricted to functions of x, sort of implying position. They're probably more 
often used in analyzing functions of time. If you're presented with a complicated function of time, how 
do you analyze it? What information is present in it? If that function of time is a sound wave you 
may choose to analyze it with your ears, and if it is music, the frequency content is just what you will 
be listening for. That's Fourier analysis. The Fourier transform of the signal tells you its frequency 
content, and sometimes subtle periodicities will show up in the transformed function even though they 
aren't apparent in the original signal. 

A function of time is f{t) and its Fourier transform is 



dtf{t)e 



iujt 



with 



fit) 



du 

27T 



g{uj)e 



-iujt 



The sign convention in these equations appear backwards from the one in Eq. (15.5), and it is. One 
convention is as good as the other, but in the physics literature you'll find this pairing more common 
because of the importance of waves. A function e*(^a;-ajf) represents a wave with (phase) velocity tu/k, 
and so moving to the right. You form a general wave by taking linear combinations of these waves, 
usually an integral. 

Example 

When you hear a musical note you will perceive it as having a particular frequency. It doesn't, and 
if the note has a very short duration it becomes hard to tell its* pitch. Only if its duration is long 



Think of a hemisemidemiquaver played at tempo prestissimo. 
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enough do you have a real chance to discern what note you're hearing. This is a reflection of the facts 
of Fourier transforms. 

If you hear what you think of as a single note, it will not last forever. It starts and it ends. Say 
it lasts from t = — T to t = +T, and in that interval it maintains the frequency Uq. 



fit) = Ae-^^"^ {-T <t< T) 
The frequency analysis comes from the Fourier transform. 



(15.11) 



T 



dte^^^ Ae-'^<'^ = A- 



T 



2A 



s'm{u — ujq)T 
{u - Uq) 



This is like the function of Eq. (15.6) except that its center is shifted. It has a peak at lu = Uq instead 
of at the origin as in that case. The width of the function g is determined by the time interval T. As T 
is large, g is narrow and high, with a sharp localization near uq. In the reverse case of a short pulse, the 
range of frequencies that constitute the note is spread over a wide range of frequencies, and you will 
find it difficult to tell by listening to it just what the main pitch is supposed to be. This figure shows 
the frequency spectrum for two notes having the same nominal pitch, but one of them lasts three times 
as long as the other before being cut off. It therefore has a narrower spread of frequencies. 

T larger 




Example 

Though you can do these integrals numerically, and when you are dealing with real data you will have 
to, it's nice to have some analytic examples to play with. I've already shown, Eq. (15.7), how the 
Fourier transform of a Gaussian is simple, so start from there. 



If g{uj) = e-(^-'^o)V'^' then 



fit) 



a 



If there are several frequencies, the result is a sum. 



In a more common circumstance you will have the time series, f{t), and will want to obtain the frequency 
decomposition, g{ui), though for this example I worked backwards. The function of time is real, but 
the transformed function g is complex. Because / is real, it follows that g satisfies g{—U!) = g*{uj). 
See problem 15.13. 
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This example has four main peaks in the frequency spectrum. The real part of g is an even 
function and the imaginary part is odd. 

^Real 





Ima! 



This is another example with four main peaks. 

In either case, if you simply look at the function of time on the left it isn't obvious what sort of 
frequencies are present. That's why there are standard, well-developed computer programs to do the 
Fourier analysis. 

15.4 Derivatives 

There are a few simple, but important relations involving differentiation. What is the Fourier transform 
of the derivative of a function? Do some partial integration. 



dt e 



iujt 



dt 



fit) 



luj / dte'^^'fit) = -iujT{f) 



(15.12) 



Here I've introduced the occasionally useful notation that F{f) is the Fourier transform of /. The 
boundary terms in the partial integration will go to zero if you assume that the function / approaches 
zero at infinity. 

The n^^ time derivative simply give you more factors: {—iu)"^ on the transformed function. 
15.5 Green's Functions 

This technique showed up in the chapter on ordinary differential equations, section 4.6, as a method to 
solve the forced harmonic oscillator. In that instance I said that you can look at a force as a succession 
of impulses, as if you're looking at the atomic level and visualizing a force as many tiny collisions by 
atoms. Here I'll get to the same sort of result as an application of transform methods. The basic 
technique is to Fourier transform everything in sight. 

The damped, forced harmonic oscillator differential equation is 



^-n^ + o-rr + kx = i'o(t) 
dt^ dt 



(15.13) 



Multiply by e*'^* and integrate over all time. You do the transforms of the derivatives by partial 
integration as in Eq. (15.12). 



dte^^^ [Eq. (15.13)] = —muj'^x — ihux + kx = Fq, where x{ijj) 
This is an algebraic equation that is easy to solve for the function x{u). 



dte^^'xit) 



x(u) 



—mu"^ — ibu + k 
Now use the inverse transform to recover the function x{t). 



x{t) 



duj 



oo 

du 



27r 



x{uj) 

iujt 



du 



-iuit 



2tt 



-mu"^ — ibu + k 



271 —rau"^ — ibu + k 



dt'Fo{t')e 



iujt 



dt'Foit'] 



du 



-itot 



271 —mu"^ — ibu + k 



(15.14) 
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In the last line I interchanged the order of integration, and in the preceding line I had to be sure to use 
another symbol t' in the second integral, not t. Now do the uj integral. 



-iujt 



27r —muj'^ — ibco + k 



-iuj{t-t') 



_^ 271 — mu;2 _ ijj^ _|_ 



(15.15) 



To do this, use contour integration. The singularities of the integrand are at the roots of the denomi- 
nator, — mcj^ — ibu + k = 0. They are 



CO 



-ib ± a/ —6^ + Akm 
2m 



Ci 

— Ti- 



( 



Both of these poles are in the lower half complex plane. The contour integral Ci is along the real axis, 
and now I have to decide where to push the contour in order to use the residue theorem. This will be 
governed by the exponential, e~*'^(*~* \ 

First take the case t < f , then e~*'^(*"* ) is of the form e+*'^, so in the complex w-plane its 
behavior in the ±i directions is as a decaying exponential toward +i (oc e"''^'). It is a rising exponential 
toward (oc e^l"^!). This means that pushing the contour Ci up toward C2 and beyond will make 
this integral go to zero. I've crossed no singularities, so that means that Eq. (15.15) is zero for t < t' . 

Next, the case that t > t' . Now e^*^(*^* ) is of the form e~*^, so its behavior is reversed from 
that of the preceding paragraph. It dies off rapidly toward —ioo and rises in the opposite direction. 
That means that I must push the contour in the opposite direction, down to C3 and to C4. Because 
of the decaying exponential, the large arc of the contour that is pushed down to —ioo gives zero for its 
integral; the two lines that parallel the i-axis cancel each other; only the two residues remain. 





C3 








du 



-iLU(t-t') 



2tt —muj"^ — ibu + k 



-27ri ^ Res 



(15.16) 



L0± 





C4 


(X) 


(X) 



The denominator in Eq. (15.15) is —m{u — io+){u — U-). Use this form to compute the residues. 
Leave the l/2n aside for the moment and you have 



-iuj{t-t') 



-iu){t-t') 



—muj'^—ibu + k —m{u — ijJ+){uj — U-) 
The residues of this at u± are the coefficients of these first order poles. 

Q—iu}+(t—t') Q—iu}-{t—t') 



at UJy. 



and at uj- 



-m[uo^ 



■ 00 A 
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The explicit values of uj± are 



2m 



and (jj. 



-ih — \/ —62 + 4A;77i 



Let 



and 



7 



2m 
6 



2m 2m 
The difference that appears in the preceding equation is then 

CJ-t- — (jj- = {(jj' — i'j) — {—u' — 27) = 2u' 

Eq. (15.16) is then 



-iuj{t-t') 



_^ 271 — ma;2 _ ^fetj + k 



,-i{uj' -i'y)(t-t') p-i{-uj'-i'y){t-t') 



-2muj' 



-l(t-t') 



+2mu' 



iuj'{t-t') _|_ 



X 



2muj' 

= -Le-^(*-*') sin U{t - t')) 
Put this back into Eq. (15.14) and you have 

(t)= / dt' Fo{t')G{t-t'), where G(t - t') = -^e''^'^*-*') sin - t')) (15.17) 

J — 00 TflUl 

If you eliminate the damping term, setting 6 = 0, this is exactly the same as Eq. (4.34). The integral 
stops att' = t because the Green's function vanishes beyond there. The motion at time t is determined 
by the force that was applied in the past, not the future. 

Example 

Apply a constant external force to a damped harmonic oscillator, starting it at time t = and keeping 
it on. What is the resulting motion? 

(t < 0) 
F, (t>0) 

where Fi is a constant. The equation (15.17) says that the solution is {t > 0) 



x{t) 



f dt' Fo{t')Git - t') = f dt' FiG{t - t') 

J -00 Jo 

sm {u'{t-t')) 



Fi / dt 



' ^_p-7(t-i'), 



2imuj' Jo 
Fi 



mcu' 

- f dt'e-^^'-''\e'^'^'-'"^ 



-iLu'{t-t')-\ 



1 



,{-7+ia;')(t-t') 



2imu' [7 — iuj' 

Fi 
2imu' 
Fi 



7 + iuo' 



^{-j~iuj'){t-t') 



t'=t 



t'=0 



2iu}' e~"f^ 

K jTT [2z7 sin u't + 2iu' cos u't] 

72 + uj'^ 72 + L ' J 

7 



m(72 + a;'2) 



1-e '^^Tcosw't + ^ sinw't 



Check the answer. If t = it is correct; x(0) = as it should. 

If t —7- 00, x{t) goes to Fi/ (m(72 + a;'2)) ; is this correct? Check it out! And maybe simplify the result 
in the process. Is the small time behavior correct? 
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15.6 Sine and Cosine Transforms 

Return to the first section of this chapter and look again at the derivation of the Fourier transform. It 
started with the Fourier series on the interval —L < x < L and used periodic boundary conditions to 
define which series to use. Then the limit as L oo led to the transform. 

What if you know the function only for positive values of its argument? If you want to write 
f{x) as a series when you know it only for < x < L, it doesn't make much sense to start the way 
I did in section 15.1. Instead, pick the boundary condition at x = carefully because this time the 
boundary won't go away in the limit that L — > oo. The two common choices to define the basis are 

m(0) = = m(L), and m'(0) = = (15.18) 

Start with the first, then Un{x) = sm{nTCx/ L) for positive n. The equation (15.2) is unchanged, save 
for the limits. 



f{x) = ^anUn{x), and {Um, f) = {Um,^anUn) = am{Um,Um) 
1 

In this basis, (um^Um) = L/2, so 



1 n=i 



oo 



(Un,f) 2 
fix) = Yl / \ Un{x) = jYl f)^n{x) 



n=l \ "/ n=l 



Now explicitly use the sine functions to finish the manipulation, and as in the work leading up to 
Eq. (15.3), denote kn = nn/ L, and the difference /^kn = vr/L. 



„ oo „r / 
2 v-^ / , , . , ,, riTix mix 
y ^ j ax j(x ) sm sm 



sin Afcn / dx'f{x')smn7rx'/L (15.19) 



TT 



For a given value of k, define the integral 

gj^{k) = / dx' sm{kx')f{x') 
Jo 

If the function / vanishes sufficiently fast as x' — > oo, this integral will have a limit as L — > oo. Call that 
limit g{k). Look back at Eq. (15.19) and you see that for large L the last factor will be approximately 
g{kn), where the approximation becomes exact as L — > oo. Rewrite that expression as 

2 °° 



fix) --Y] sm{knX)Akn gikn) (15.20) 

TT ^ 

As L — > OO, you have Afc^ — > 0, and that turns Eq. (15.20) into an integral. 

2 /"OO /"OO 

f(x) = — / dks'mkxgik), where gik) = I dxsmkxf(x) (15.21) 
Jo Jo 
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This is the Fourier Sine transform. For a parallel calculation leading to the Cosine transform, see 
problem 15.22, where you will find that the equations are the same except for changing sine to cosine. 

f(x) = — / dkcoskxgik), where 9{k) = / dxcoskxf(x) (15.22) 
71" Jo Jo 

What is the sine transform of a derivative? Integrate by parts, remembering that / has to 
approach zero at infinity for any of this to make sense. 

POO oo /.oo POO 

/ dx sinkx f'{x) = smkxf{x) ~ k dx coskx f{x) = —k / dx coskx f{x) 
Jo Jo Jo 

For the second derivative, repeat the process. 

roo roo 

/ dx sinkx f"{x) = kf{0) - k"^ dx sinkx f{x) (15.23) 
Jo Jo 



15.7 Wiener-Khinchine Theorem 

If a function of time represents the pressure amplitude of a sound wave or the electric field of an 
electromagnetic wave the power received is proportional to the amplitude squared. By Parseval's identity, 
the absolute square of the Fourier transform has an integral proportional to the integral of this power. 
This leads to the interpretation of the transform squared as some sort of power density in frequency. 
\g{u)\'^du is then a power received in this frequency interval. When this energy interpretation isn't 
appropriate, |(7(t<;)p is called the "spectral density." A useful result appears by looking at the Fourier 
transform of this function. 



du 
27r 



2„—iuit 



du 
2tx 



-iu)t 



g-[uj)e I dt'f{t')e 



i\Au)t' 



duj 



dt'fit') ^9*{uj)e 



Aojt' „—iujt 



dt'fit') 



27T 

du 



271 

dt'fit')f{t'-t) 



^(u;)e-^^(^'-*) 



(15.24) 



When you're dealing with a real /, this last integral is called the autocorrelation function. It tells you 
in some average way how closely related a signal is to the same signal at some other time. If the signal 
that you are examining is just noise then what happens now will be unrelated to what happened a few 
milliseconds ago and this autocorrelation function will be close to zero. If there is structure in the signal 
then this function gives a lot of information about it. 

The left side of this whole equation involves two Fourier transforms ( / — > then to it's 
transform). The right side of this theorem seems to be easier and more direct to compute than the left, 
so why is this relation useful? It is because of the existence of the FFT, the "Fast Fourier Transform," 
an algorithm that makes the process of Fourier transforming a set of data far more efficient than doing 
it by straight-forward numerical integration methods — faster by factors that reach into the thousands 
for large data sets. 
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Problems 

15.1 Invert the Fourier transform, g, in Eq. (15.7). 

15.2 What is the Fourier transform of e^k^x-x'^ /o'^-f /^^^. translation of the /cq = case 

2 / 2 

15.3 What is the Fourier transform of a;e~^ ? 

15.4 What is the square of the Fourier transform operator? That is, what is the Fourier transform of 
the Fourier transform? 



15.5 A function is defined to be 

fix) 



1 {—a < X < a) 
(elsewhere) 



What is the convolution of / with itself? (/ * /)(x) And graph it of course. Start by graphing both 
fix') and the other factor that goes into the convolution integral. 
Ans: (2a — |x|) for (—2a < x < +2a), and zero elsewhere. 

15.6 Two functions are 

„.._fl {a<x <h) f/\/l {A<x < B) 

JiW-|q (elsewhere) /2W - | q (elsewhere) 

What is the convolution of /i with /2? And graph it. 

15.7 Derive these properties of the convolution: 

{a) f*g = g*f {h) f*{g*h) = if*g)*h {c)6{f*g) = f *5g+g*6f where 6 f{t) =tf{t), 
Sg{t) = tg{t), etc. (d) What are 6^{f * g) and 5^{f * g)l 

15.8 Show that you can rewrite Eq. (15.9) as 

Hf*g)=Tif)-Tig) 

where the shorthand notation T{f) is the Fourier transform of /. 

15.9 Derive Eq. (15.10) from Eq. (15.9). 

15.10 What is the analog of Eq. (15.10) for two different functions? That is, relate the scalar product 
of two functions. 



(/l,/2>= / /;(X)/2(X) 

J —00 

to their Fourier transforms. Ans: J gl{k)g2ik) dk/2TT 

15.11 In the derivation of the harmonic oscillator Green's function, and starting with Eq. (15.15), I 
assumed that the oscillator is underdamped: that 6^ < Akm. Now assume the reverse, the overdamped 
case, and repeat the calculation. 

15.12 Repeat the preceding problem, but now do the critically damped case, for which = Akm. 
Compare your result to the result that you get by taking the limit of critical damping in the preceding 
problem and in Eq. (15.17). 
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15.13 Show that if f{t) is real then the Fourier transform satisfies g{—ui) = g*{uj). 
What are the properties of if / is respectively even or odd? 



15.14 Evaluate the Fourier transform of 

fix)- 



A{a — \x\) {—a<x<a) 
(otherwise) 

How do the properties of the transform vary as the parameter a varies? 
Ans: 2A{1 - coska)/P 

15.15 Evaluate the Fourier transform of Ae~^^^^. Invert the transform to verify that it takes you back 
to the original function. Ans: 2a/ {a^ + k"^) 

15.16 Given that the Fourier transform of f{x) is g{k), what is the Fourier transform of the function 
translated a distance a to the right, fi{x) = f{x — a)l 



15.17 Schroedinger's equation is 



dt ~ 2m dx^ + ^ ^^^^ 



Fourier transform the whole equation with respect to x, and find the equation for ^{k,t), the Fourier 
transform 'ip{x,t). The result will not be a differential equation. 
Ans: -ihd<^{k, t) /dt = {h'^P/2m)<^ + (v * $) /2tt 



15.18 Take the Green's function solution to Eq. (15.13) as found in Eq. (15.17) and take the limit 
as both k and b go to zero. Verify that the resulting single integral satisfies the original second order 
differential equation. 

15.19 (a) In problem 15.18 you have the result that a double integral (undoing two derivatives) can 
be written as a single integral. Now solve the equation 



d^x 



Fit) 



C 



directly, using the same method as for Eq. (15.13). You will get a pole at the origin and how do 
you handle this, where the contour of integration goes straight through the origin? Answer: Push the 
contour up as in the figure. Why? This is what's called the "retarded solution" for which the value 
of x{t) depends on only those values of F{t') in the past. If you try any other contour to define the 
integral you will not get this property. (And sometimes there's a reason to make another choice.) 
(b) Pick a fairly simple F and verify that this gives the right answer. 
Ans: lf_^dfF{t'){t-t')^ 

15.20 Repeat the preceding problem for the fourth derivative. Would you care to conjecture what 31/2 
integrals might be? Then perhaps an arbitrary non-integer order? 
Ans: lf_^dt'F{t'){t-tr 



15.21 What is the Fourier transform of xf{x)l Ans: ig'ik) 
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15.22 Repeat the calculations leading to Eq. (15.21), but for the boundary conditions u'{0) = = 
u'{L), leading to the Fourier cosine transform. 

15.23 For both the sine and cosine transforms, the original function f{x) was defined for positive x 
only. Each of these transforms define an extension of / to negative x. This happens because you 
compute g{k) and from it get an inverse transform. Nothing stops you from putting a negative value 
of X into the answer. What are the results? 

15.24 What are the sine and cosine transforms of e~'^^ . In each case evaluate the inverse transform. 

15.25 What is the sine transform of f{x) = 1 for < x < L and f{x) = otherwise. Evaluate the 
inverse transform. 

15.26 Repeat the preceding calculation for the cosine transform. Graph the two transforms and com- 
pare them, including their dependence on L. 

15.27 Choose any different way around the pole in problem 15.19, and compute the difference between 
the result with your new contour and the result with the old one. Note: Plan ahead before you start 
computing. 



